First results from the CERN Axion Solar Telescope (CAST) 



K. Zioutas,® S. Andriamonje,^ V. Arsov,^^-"* S. Aune,^ D. Autiero,^'0 F. Avignone,^ K. Barth,^ A. Belov," 
B. Beltraii,^ H. Brauninger,^ J. M. Carmona,^ S. Cebriaii,^ E. Chesi,i J. I. Collar/ R. Creswick,^ T. Dafni,^ 
M. Davenport,^ L. Di Leila, -^'0 C. Eleftheriadis,® J. Englhauser,^ G. Fanourakis,^ H. Farach,'^ E. Ferrer,^ 
H. Fischer, J. Franz, P. Friedrich,^ T. Geralis,^ I. Giomataris,^ S. Gninenko,^^ N. Goloubev,^^ 
M. D. Hasinoff,^^ F. H. Heinsius,^" D.H.H. Hoffmann,"* I. G. Irastorza,^ J. Jacoby,*^ D. Kang,^° K. Konigsniann. 
R. Kotthaus,*'* M. Krcmar,*^ K. Kousouris,^ M. Kuster,^ B. Lakic,^^ C. Lasseur,* A. Liolios,^ A. Ljubicic,*^ 
G. Lutz,!-* G. Luzon,^ D. W. Miller,'^ A. Morales,*^'|l J. Morales,^ M. Mutterer," A. Nikolaidis,^ A. Ortiz,^ 
T. Papaevangelou,! A. Placet ^ G. Raffelt,!'^ J. Ruz,*^ H. Riege,^ M. L. Sarsa,^ I. Savvidis,^ W. Serber,** 
P. Serpico,*'* Y. Semertzidis,^^|j L. Stewart,^ J. D. Vieira,^ J. Villar,^ L. Walckiers,* and K. Zachariadou^ 

(CAST Collaboration) 
^European Organization for Nuclear Research (CERN), Geneve, Switzerland 
^ DAPNIA, Centre d'Etudes Nucleaires de Saclay (CEA-Saclay), Cif-sur-Yvette, France 
^Department of Physics and Astronomy, University of South Carolina, Columbia, SC, USA 
CSI-Darmstadt and Institut fiir Kernphysik, TU Darmstadt, Darmstadt, Germany 
Max- Planck- Institut fiir Extraterrestrische Physik, G arching, Germany 
^Institute de Fisica Nuclear y Altas Energias, Universidad de Zaragoza, Zaragoza, Spain 
''Enrico Fermi Institute and KICP, University of Chicago, Chicago, IL, USA 
^Aristotle University of Thessaloniki, Thessaloniki, Greece 
''National Center for Scientific Research "Demokritos" , Athens, Greece 
^" Albert-Ludwigs-Universitdt Freiburg, Freiburg, Germany 
^'institute for Nuclear Research (INR), Russian Academy of Sciences, Moscow, Russia 
'^Department of Physics and Astronomy, University of British Columbia, Vancouver, Canada 
Johann Wolfgang Goethe- Universitdt, Institut fiir Angewandte Physik, Frankfurt am Main, Germany 
'^^^ Max- Planck- Institut fiir Physik (Wemer-Heisenberg-Institut) , Munich, Germany 
'^Rudjer Boskovic Institute, Zagreb, Croatia 
(Dated: February 7, 2008) 

Hypothetical axion-like particles with a two-photon interaction would be produced in the Sun 
by the Primakoff process. In a laboratory magnetic field ( "axion helioscope" ) they would be trans- 
formed into X-rays with energies of a few keV. Using a decommissioned LHC test magnet, CAST ran 
for about 6 months during 2003. The first results from the analysis of these data are presented here. 
No signal above background was observed, implying an upper limit to the axion-photon coupling 
ga-i < 1.16 X 10~^° GeV~^ at 95% CL for rua 0.02 eV. This limit, assumption-free, is comparable 
to the limit from stellar energy-loss arguments and considerably more restrictive than any previous 
experiment over a broad range of axion masses. 

PACS numbers: 95.35.-|-d; 14.80.Mz; 07.85. Nc; 84.71. Ba 
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Introduction. — Neutral pions, gravitons, hypothetical 
axions, or other particles with a two-photon interac- 
tion can transform into photons in external electric or 
magnetic fields, an effect first discussed by Primakoff 
in the early days of pion physics 0]. Therefore, stars 
could produce these particles by transforming thermal 
photons in the fluctuating electromagnetic fields of the 
stellar plasma In laboratory or astrophysical 

B-fields, transitions between these particles and pho- 
tons occur 0, 3, an effect that can be observed in 
the laboratory a|, affects the propagation of cosmic 7- 
rays 0, @], and can modify the app arent brightness of 
distant astronomical sources HEJini- Gravitons in- 
teract too weakly to be observable in these situations, 
while pions are too heavy. However, these effects can be 
crucial for new particles, notably the pseudoscalar ax- 
ions that arise in the context of the Peccei-Quinn so- 
lution to the strong CP problem and are viable cold 



dark matter candidates 0, ^| . Galactic axions are cur- 
rently being sought by two large-scale Primakoff-type 
microwave cavity experiments |l2l |. Anomalous stel- 
lar energy loss by axion emission is constrained by the 
observed properties of globular cluster stars, implying 
Qa-y ^ 10~^° GeV~^ [3 for the axion-photon coupling, 
where the axion-photon interaction is written in the usual 
form Ca-y = — jga-yF^vF^'^ a = ga-y E • B a. Axions would 
also contribute to the magnetically induced vacuum bire- 
fringence, interfering with the corresponding QED ef- 
fect jlllil. The PVLAS experiment [l3| apparently ob- 
serves such an effect far in excess of the QED expectation, 
although an interpretation in terms of axion-like particles 
requires a coupling strength far larger than existing lim- 
its. 

On the other hand, the Sun would be a strong ax- 
ion source and thus offers a unique opportunity to ac- 
tually detect these particles by taking advantage of 
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their back-conversion into X-rays in laboratory mag- 
netic fields 0. The expected solar axion flux at 
the Earth due to the Primakoff process is — 
glo 3.67 X 10" cm'^ s"! (where 510 = gaj GeV) 
with an approximate spectrum d^a/dEa — gfg 3.821 x 
IQio cm-2 s-i keV-i (£;,/keV) V(e^"/i i03 kcV _ 
an average energy of 4.2 keV [l5j. Axion interactions 
other than the two-photon vertex would provide for ad- 
ditional production channels, but in the most interesting 
scenarios these channels are severely constrained, leaving 
the Primakoff effect as the dominant one In any case, 
it is conservative to use the Primakoff effect alone when 
deriving limits on ga-y. 

The conversion probability in a -B-field in vacuum is 
Pa^-y — (gajB/q)'^ sm^{qL/2), where L is the path 
length and q = m'^/2Ea is the axion-photon momentum 
difference. For qL ^ 1 where the axion-photon oscilla- 
tion length far exceeds L we have Pa^j — {gayBL/2)'^, 
implying an X-ray flux of 

*, = O.Mc„-d-.A(^)'(^)'. (1, 

For qL ^ 1 this rate is reduced due to the axion- 
photon momentum mismatch. The presence of a gas 
would provide a refractive photon mass so that 
q = — ■m^\/2Ea. For nia ~ the maximum rate 
can thus be restored . 

The first implementation of the axion helioscope con- 
cept was performed in More recently, the Tokyo 
axion helioscope of i = 2.3 m and i? = 3.9 T has 
provided the Hmit 510 < 6.0 at 95% CL for rua < 0.03 eV 
(vacuum) and gio < 6.8-10.9 for < 0.3 eV (using a 
variable- pressure buffer gas) [i^ . Limits from crystal de- 
tectors 113, 13 ^-^s much less restrictive. 

CAST Experiment. — In order to detect solar axions 
or to improve the existing limits on gaj an axion he- 
lioscope has been built at CERN by refurbishing a de- 
commissioned LHC test magnet |25j which produces a 
magnetic field oi B = 9.0 T in the interior of two par- 
allel pipes of length L ~ 9.26 m and a cross-sectional 
area yl = 2 x 14.5 cm^. The aperture of each of the 
bores fully covers the potentially axion-emitting solar 
core 1/lOth of the solar radius). The magnet is 
mounted on a platform with ±8° vertical movement, al- 
lowing for observation of the Sun for 1.5 h at both sunrise 
and sunset. The horizontal range of ±40° encompasses 
nearly the full azimuthal movement of the Sun through- 
out the year. The time the Sun is not reachable is devoted 
to background measurements. A full cryogenic station is 
used to cool the superconducting magnet down to 1.8 K 
[iif. The hardware and software of the tracking system 
have been precisely calibrated, by means of geometric 
survey measurements, in order to orient the magnet to 
any given celestial coordinates. The overall CAST point- 
ing precision is better than 0.01° 0| . At both ends of the 



magnet, three different detectors have searched for excess 
X-rays from axion conversion in the magnet when it was 
pointing to the Sun. Covering both bores of one of the 
magnet's ends, a conventional Time Projection Chamber 
(TPC) is looking for X-rays from "sunset" axions. At 
the other end, facing "sunrise" axions, a second smaller 
gaseous chamber with novel MICROMEGAS (micromesh 
gaseous structure - MM) readout is placed behind 
one of the magnet bores, while in the other one a focus- 
ing X-ray mirror telescope is working with a Charge Cou- 
pled Device (CCD) as the focal plane detector. Both the 
CCD and the X-ray telescope are prototypes developed 
for X-ray astronomy jl^. The X-ray mirror telescope 
can produce an "axion image" of the Sun by focusing 
the photons from axion conversion to a 6 mm^ spot on 
the CCD. The enhanced signal-to-background ratio sub- 
stantially improves the sensitivity of the experiment. A 
detailed account of the technical aspects of the experi- 
ment will be given elsewhere. 

Data Analysis and First Results. — CAST operated for 
about 6 months from May to November in 2003, during 
most of which time at least one detector was taking data. 
The results presented here were obtained after the anal- 
ysis of the data sets listed in Tabled An independent 
analysis was performed for each data set. Finally, the 
results from all data sets are combined. 

An important feature of the CAST data treatment is 
that the detector backgrounds are measured with '^10 
times longer exposure during the non-alignment periods. 
The use of these data to estimate and subtract the true 
experimental background during Sun tracking data is the 
most sensitive step in the CAST analysis. To assure the 
absence of systematic effects, the main strategy of CAST 
is the use of three independent detectors with comple- 
mentary approaches. In the event of a positive signal, 
it should appear consistently, if strong enough, in each 
of the three detectors when it is pointing at the Sun. 
In addition, an exhaustive recording of experimental pa- 
rameters was done, and a search for possible background 
dependencies on these parameters was performed. A de- 
pendence of the TPC background on the magnet po- 
sition was found, caused by its relatively large spatial 
movements at the far end of the magnet, which resulted 
in appreciably different environmental radioactivity lev- 
els. Within statistics, no such effect was observed for the 
sunrise detectors which undergo a much more restricted 
movement. To correct for this systematic effect in the 
TPC data analysis, an effective background spectrum is 
constructed only from the background data taken in mag- 
net positions where Sun tracking has been performed and 
this is weighted accordingly with the relative exposure of 
the tracking data. Further checks have been performed 
in order to exclude any possible systematic effect. They 
were based on rebinning the data, varying the fitting win- 
dow, splitting the data into subsets and verifying the null 
hypothesis test in energy windows or areas of the detec- 
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TABLE L Data sets included in our result. 



Data set 


Tracking exposure Background exposure {s*^)bcstflt (ilf error) 


xLll/d-o.f 


■ /d.o.f 






(h) 


(h) 


(10-« GeV-4) 






(10-10 GeV-i) 


TPC 


62.7 


719.9 


-1.1 ± 3.3 


18.2/18 


18.1/17 


1.55 


MM set A 


43.8 


431.4 


-1.4 ±4.5 


12.5/14 


12.4/13 


1.67 


MM set B 


11.5 


121.0 


2.5 ±8.8 


6.2/14 


6.1/13 


2.09 


MM set C 


21.8 


251.0 


-9.4 ±6.5 


12.8/14 


10.7/13 


1.67 


CCD 


121.3 


1233.5 


0.4 ±1.0 


28.6/20 


28.5/19 


1.23 



xlO 
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FIG. 1: Panels (a) and (b) show respectively the experimental 
subtracted spectrum of the TPC data set and MM data set 
A, together with the expectation for the best fit Qa-y (lower 
curve) and for the 95% CL limit on gaj- For (a) the verti- 
cal dashed lines indicate the fitting window. The structure 
at 3 keV in the expected spectrum reflects the change in the 
efficiency curves due to the Ar K-edge of the detector gas 
mixtures. Panel (c) shows both the tracking (dots) and back- 
ground (dashed line) spectra of the CCD data set, together 
with the expectation (background plus signal) for ga-y at its 
95% CL limit, in units of total counts in the restricted CCD 
area (54.3 mm'^) in the tracking exposure time (121.3 h). 



tors where no signal is expected. In general, the system- 
atic uncertainties are estimated to have an effect of less 
than ^10% in the final upper limits obtained. 

For a fixed ma, the theoretically expected spectrum 



of axion-induced photons has been calculated and multi- 
plied by the detector efficiency curves (determined both 
by calculation and measurements). These spectra, which 
are proportional to g^^, are directly used as fit functions 
to the experimental subtracted spectra (tracking minus 
background) for the TPC and MM. For these data, the 
fitting is performed by standard minimization. Re- 
garding the CCD data, the analysis is restricted to the 
small area on the CCD where the axion signal is expected 
after the focusing of the X-ray telescope. During the 
data taking period of 2003 a continuous monitoring of 
the pointing stability of the X-ray telescope was not yet 
possible, therefore a signal area larger than the size of the 
sun spot had to be considered. Taking into account all 
uncertainties of the telescope alignment, the size of the 
area containing the signal was conservatively estimated 
to be 34 X 71 pixels (54.3 mm^). As in the other detec- 
tors, the background is defined by the data taken from 
the same area during the non-tracking periods but, al- 
ternatively, the background in the signal area was also 
determined by extrapolating the background measured 
during tracking periods in the part of the CCD not con- 
taining the sun spot. Both methods of background se- 
lection led to the same final upper limit on the coupling 
constant ga-y- The resulting low counting statistics in 
the CCD required the use of a likelihood function in the 
minimization procedure, rather than a x^-analysis. 

The best fit values of obtained for each of the data 
sets are shown in Table together with their la error 
and the corresponding Xmin values and degrees of free- 
dom. Figure ^ shows the plots of some of those fits. 
Each of the data sets is individually compatible with the 
absence of any signal as can be seen from the XnuU values 
shown in Table HJ The excluded value of g^^ was con- 
servatively calculated by taking the limit encompassing 
95% of the physically allowed part (i.e. positive signals) 
of the Bayesian probability distribution with a flat prior 
in g^-^. The described procedures were done using g^^ in- 
stead of ga-y as the minimization and integration param- 
eter because the signal strength (i.e. number of counts) 
is proportional to g^^. The 95% CL limits on gaj for 
each of the data sets are shown in the last column of 
Table m They can be statistically combined by multiply- 
ing the Bayesian probability functions and repeating the 
previous process to find the combined result for the 2003 
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CAST data: 



ga-y < 1.16 X 10~^°GeV~^(95%CL). 

Thus far our analysis was limited to the mass range 
TTia ^ 0.02 eV where the expected signal is mass- 
independent because the axion-photon oscillation length 
far exceeds the length of the magnet. For higher iria the 
overall signal strength diminishes rapidly and the spec- 
tral shape differs. Our procedure was repeated for differ- 
ent values of to obtain the entire 95% CL exclusion 
line shown in Fig. |21 

Summary. — Our limit improves the best previous lab- 
oratory constraints |2C| on ga^y by a factor 5 in our co- 
herence region nia ^ 0.02 eV. This result excludes an 
important part of the parameter space not excluded by 
solar age considerations l30l| and is comparable, in this 
range of masses, to the limit derived from stellar energy- 
loss arguments. A higher sensitivity is expected from 
the 2004 data with improved conditions in all detectors, 
which should allow us to surpass the astrophysical limit. 
In addition, starting in 2005, CAST plans to take data 
with a varying-pressure buffer gas in the magnet pipes, 
in order to restore coherence for axion masses above 0.02 
eV. The extended sensitivity to higher axion masses will 
allow us to enter into the region shown in Fig. [5] which is 
especially motivated by axion models j3ll |. 

Acknowledgments. — We thank CERN for hosting the 
experiment and for the contributions of J. P. Bojon, 
F. Cataneo, R. Campagnolo, G. CipoUa, F. Chiusano, 
M. Delattre, F. Formenti, M. Genet, J. N. Joux, A. Lip- 
pitsch, L. Musa, R. De Oliveira, A. Onnela, C. Ros- 
set and B. VuUierme. Wc thank in particular F. James 
for his advice concerning the statistical treatment of the 
data. We acknowledge support from NSERC (Canada), 
MSES (Croatia), CEA (France), BMBF (Germany), 
GSRT (Greece), RFFR (Russia), CICyT (Spain), and 
NSF (USA). We dedicate this paper to the memory of 
our friend and colleague Angel Morales. 



* Present addr.: Inst, de Physique Nucleaire, Lyon, France. 
^ Present addr.; Scuola Normale Superiore, Pisa, Italy. 

* Deceased. 

^ Permanent addr.; Brookhaven Nat. Lab., NY-USA. 
[1] H. Primakoff, Phys. Rev. 81, 899 (1951). 
[2] D. A. Dicus et ai, Phys. Rev. D 18, 1829 (1978). 
[3] G. G. Raffelt, Ann. Rev. Nucl. Part. Sci. 49, 163 (1999). 
[4] P. Sikivie, Phys. Rev. Lett. 51, 1415 (1983) [Erratum- 

ibid. 52, 695 (1984)]. 
[5] G. Raffelt and L. Stodolsky, Phys. Rev. D 37, 1237 

(1988). 

[6] S. Eidelman et al, Phys. Lett. B 592, 1 (2004). 

[7] D. S. Gorbunov, G. G. Raffelt and D. V. Semikoz, Phys. 

Rev. D 64, 096005 (2001). 
[8] C. Csaki et al, JCAP 0305, 005 (2003). 




w 



10' 10"^ 10'^ 10'^ 10'^ 1 10 

maxion(eV) 

FIG. 2; Exclusion limit (95% CL) from the CAST 2003 data 
compared with other constraints discussed in the introduc- 
tion. The shaded band represents typical theoretical models. 
Also shown is the future CAST sensitivity as foreseen in the 
experiment proposal. 



[9] C. Csaki, N. Kaloper and J. Terning, Phys. Rev. Lett. 
88, 161302 (2002). 

[10] E. Mortsell and A. Goobar, JCAP 0304, 003 (2003). 

[11] B. A. Bassett, Astrophys. J. 607, 661 (2004). 

[12] R. Bradley et al, Rev. Mod. Phys. 75, 777 (2003). 

[13] L. Maiani, R. Petronzio and E. Zavattini, Phys. Lett. B 
175, 359 (1986). 

[14] G. Cantatore et al, presented at IDM 2004, Edinburgh, 
England, 4-10 Sep^004. 

[15] Th e sp ectrum in |la| has been changed to that proposed 
in however with a modified normalization constant 
to match the total axion flux used here, which is predicted 
by a more recent solar model fl^ . 

[16] R. J. Creswick et al, Phys. Lett. B 427, 235 (1998). 

[17] J. N. Bahcall, M. H. Pinsonneault and S. Basu, Astro- 
phys. J. 555, 990 (2001). 

[18] K. van Bibber et al, Phys. Rev. D 39, 2089 (1989). 

[19] D. M. Lazarus et al, Phys. Rev. Lett. 69, 2333 (1992). 

[20] S. Moriyama et al, Phys. Lett. B 434, 147 (1998). 

[21] Y. Inoue et al, Phys. Lett. B 536, 18 (2002). 

[22] F. T. Avignone et al, Phys. Rev. Lett. 81, 5068 (1998). 

[23] A. Morales et al, Astropart. Phys. 16, 325 (2002). 

[24] R. Bernabei et al, Phys. Lett. B 515 6 (2001). 

[25] K. Zioutas et al, Nucl. lustrum. Meth. A 425 480 (1999). 

[26] K. Barth et al, Proc. 2003 Cryogenic Engineering 
Conference (CEC) and Cryogenic Materials Conference 
(ICMC). 

[27] http: / /cast .web.cern.ch / CAST /edited_tracking.movl 
[28] Y. Giomataris et al, Nucl. lustrum. Meth. A 376 29 
(1996). 

[29] J. Altmann et al, in Proceedings of SPIE: X-Ray Optics, 
Instruments, and Mission, 1998, edited by Richard B. 
Hoover and Arthur B. Walker, p. 350; J. W. Egle et al, 
ibid., p. 359; P. Friedrich et al, ibid., p. 369. 

[30] H. Schlattl, A. Weiss and G. Raffelt, Astropart. Phys. 10 



353 (1999). 

D. B. Kaplan, Nucl. Phys. B 260, 215 (1985). 



